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ABSTRACT 

Using the VLA, a series of high velocity resolution observations have been made of the 
M82 starburst at 1.6 GHz. These observations follow up on previous studies of the main 
line OH maser emission in the central kiloparsec of this starburst region, but with far 
greater velocity resolution, showing significant velocity structure in some of the maser 
spots for the first time. A total of thirteen masers were detected, including all but one 
of the previously known sources. While some of these masers are still unresolved in 
velocity, these new results clearly show velocity structure in spectra from several of 
the maser regions. Position- velocity plots show good agreement with the distribution 
of Hi including interesting velocity structure on the blue-ward feature in the west of 
the starburst which traces the velocity distribution seen in the ionised gas. 
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1 INTRODUCTION 

Masers are common in star forming regions within the Milky 
Way and many nearby galaxies are now known to con- 
tain masers of one flavour or another, with powerful mega- 
masers detected out to greater and greater redshifts (e.g. 
the recent discovery of a water maser at a redshift of 2.639; 
llmpellizzeri et aT1l2008h . Maser sources require dense molec- 
ular gas and a nearby source of energy to create the required 
population inversion. The presence of one or more type of 
maser in a given region can therefore be used to probe the 
physical conditions of the gas, making them useful tools for 
understanding complex starburst systems. 

Galactic OH masers located in star forming regions are 
typically very compact (~10 13 cm), have very narrow line 
widths (~lkms _1 ) and flux ratios (1667/1665) less than 
one. In contrast, megamaser emission in ULIRGs tends to 
be spatially unresolved, have line widths >100kms _ , flux 
ratios >1 and total intensities typi cally a f actor of 10 6 
greater that that of Galactic sources (|Loll2005l ) . In between 
these classes lie the so-called "kilomasers" with luminosi- 
ties around 10 3 timer greater than typical Galactic sources. 
These are found in nearby starbursts where the star forma- 
tion rate is intermediate between that of the Milky Way and 
those of distant ULIRGs. 

The nearby starburst M82 is the prototypical starburst 
and one of the most intensively stud ied. It is known to con- 
tain large amoun ts of gas (e.g. Hi, IWills et al. 2000; CO, 
IShen fc Lolll995h , have strong outflows (e.g. winds traced 



by Ha, IShopbell fc Bland-HawthoTnl 1 19981). numerous su- 



pernova remnants and Hn regions (e.g. Muxlow et al.lll994l : 



iMcDonald et al.ll2002l ; iFenech et al.ll2008l ) and a history of 
strong star formation due to a past encounter with M81, al- 
though there is evidence that the rate o f star formation ha s 
decreased significantly in the last 5 Myr (|Beirao et alfeooih . 

Various types of maser have been detected in the cen- 
tral starburst of M82 over the last few decades, as would 
be expected for a galaxy with a significant molecular gas 
content and strong continuum emission, and these features 
can be used to investigate the molecular gas conditions and 
dynamics within th e galaxy. Main line OH masers were first 
detected in M82 bv lWeliachew et alj (|l984r i and, while they 
are significantly stronger than maser emission seen in the 
Milky Way, they are fainter than the megamasers seen at 
the other extreme in AGN. 

Previous studies of the OH masers in M82 l|Argo et al.l 

12007ft have found 10 individual masing regions, none of 
which are spatially resolved with the l'.'4 beam of the VLA. 
These observations, however, were designed to probe the OH 
absorption in M82, rather than the masers, so the spectral 
resolution was comparatively poor. 

This paper describes follow-up observations carried out 
with the VLA in 2006 at higher spectral resolution in order 
to investigate the maser population in more detail. Unless 
otherwise stated, all positions are given in J2000 coordinates 
in the format bKbb corresponding to 09 h 55 m 

and +69°40'bb'.'bb respectively. At the distance of M82 
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ObsID 


Date 


Time 


(LST) 


AA302a 


Feb 7/8 


0200 


- 1030 


AA302b 


Feb 16/17 


0230 


- 1030 


AA302c 


Feb 21 


0300 


- 0900 


AA302d 


Mar 19 


0330 


- 0930 


AA302e 


Mar 27/28 


0430 


- 1230 



Table 1. VLA observations carried out in 2006 during pro- 
gramme AA302. The structure of each observation was similar 
and included scans on flux and bandpass calibrators as well as 
scans on M82 and 0954+745 at two frequencies designed to cover 
the required bandwidth. 

f3.6Mpc, lFreedman et allll994h . one arcsecond corresponds 
to a linear distance of 17.5 parsecs. 



2 OBSERVATIONS 

High velocity resolution observations were carried out us- 
ing the VLA in A configuration during February and March 
2006. These observations were intended to to enable the de- 
tection of masers with narrow velocity widths which were 
not bright enough to be detectable in the much wider chan- 
nels of the 2002 observations, and attempt to resolve velocity 
structure in those already known. 

The experiment was carried out in five separate runs, 
details of which are given in Table[T] The sources 1331+305, 
0319+415 and 0954+745 were used as flux, bandpass and 
phase calibrators respectively. The observations each con- 
sisted of scans on the flux and bandpass calibrators as well 
as scans on M82 and 0954+745 at two frequencies such that 
the total frequency coverage was sufficient to cover the en- 
tire velocity range expected across M82 at both 1665 and 
1667 MHz with some overlap between IFs at the band edges. 
To achieve this, the correlator was used in 4-IF mode with 
a bandwidth of 1.5 MHz and 128 channels per IF, resulting 
in a velocity resolution of 2 kms -1 . Velocities in this paper 
are quoted relative to 225 kms -1 , the systemic radio LSR 
velocity of M82. 

The data from each run were edited and calibrated sep- 
arately before being combined prior to imaging. For each 
run, the process was as described below. Firstly, the channel 
zero file (created by averaging the inner 75 per cent of the 
band) was examined for each source, polarization, IF and 
frequency in order to detect major problems such as dead 
antennas and bad scans. This pseudo-continuum dataset was 
then calibrated using standard methods for continuum cali- 
bration with the VLA. Each of the four IFs were calibrated 
separately. Together with the flag table, the final calibra- 
tion tables from this step, one for each IF, were then copied 
over to the corresponding line datasets. A bandpass calibra- 
tion was then performed for each frequency using 0319+415 
before the data were split out separately for each IF. The 
datasets for each of the epochs were then combined to cre- 
ate one dataset per IF from which image cubes were made. 
Each of the four cubes were then inspected separately using 
various tasks within AIPS. 

When the first run, AA302a, was initially examined in 
the AIPS task POSSM, a large amount of interference was dis- 
covered. The data were examined more closely in the visual 
task SPFLG and strong interference was found to be present 



on all sources, varying both in time and with baseline. It is 
not telescope specific, nor is it apparent on all baselines at 
a given time. It is generally stronger at 1665 MHz than in 
the 1667 MHz data and both stronger and more prevalent in 
Stokes RR than LL, but largely it does only affect a small 
subset of the channels. A patient examination of each source 
in each frequency, stokes and IF was made in SPFLG for each 
run and data were flagged on all sources, polarisations, IFs 
and times where the problem was seen to appear. Several 
checks were made throughout this process to ensure that all 
interference was removed. 

The resulting calibrated cubes, one for each of the four 
IFs, were imaged and then continuum subtracted using a set 
of line free channels at either end of the band before being 
examined for maser features. The mean 3<r noise in each cube 
was 2.7 mjy beam -1 channel -1 (similar sensitivity to the low 
velocity resolution MERLIN data from 1995 but a factor of 
almost six worse than the 2002 VLA observations due to 
the much narrower channels). By comparing measurements 
of continuum sources in M82 it was found that, to within 
the uncertainties, the positions and flux scales matched in 
all four final datasets. 



3 RESULTS AND GENERAL COMMENTS 

Table [2] lists all the maser features in this dataset which 
were present at 5a or more in more than one channel at 
either one or both frequencies, together with their positions, 
peak fluxes or limits in each line, flux ratio between the 
two lines and the velocity at the peak of the line profile. 
Figure [T] shows the maser locations across the central kpc 
of the M82 starburst. This figure was constructed from the 
four final continuum-subtracted image cubes: each cube was 
collapsed in frequency, ignoring the five channels at either 
end of the band in each cube, and the resulting images were 
then summed to create one image. The diagram shows all of 
the maser features above 5cr in more than one consecutive 
channel, whether they are detected at 1665 MHz, 1667 MHz 
or both. 

Despite the lower sensitivity per channel than in the 
VLA 2002 dataset, all but one of the previously known 
masers appear in this dataset. Due to the narrower chan- 
nels used in these observations, the masers generally appear 
much brighter here than in the 2002 dataset. The exception 
to this is 50.02+45.8, a maser which has remained unde- 
tected since MERLIN observations in 1995 (see Section [S]). 

In general, fluxes measured in this dataset are more 
than a factor of two brighter than those measured in the 
much broader channels used in the 2002 observations. The 
smallest difference is 50.95+45.4 where the difference is a 
factor of 1.9 at 1667 MHz (2.78 at 1665 MHz), and the largest 
is 53.11+47.9 which is brighter at 1665 MHz by a factor of 
>10 in 2006. This feature is narrow even at a velocity reso- 
lution of 2 kms -1 so could be brighter still. 

Positions measured in the 2006 data are all within half 
a beam of those measured in 2002 at the same spatial reso- 
lution (l'.'4). The mean positional offset between sources in 
the two epochs is 0'.'17 with 51.94+48.3 having the largest 
offset (0'.'36). 

Spectra from each maser feature are shown in Figure 
[2] (1665 MHz) and Figure [3] (1667 MHz). Most of the maser 
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J2000 ID 


R.A. 


Dec. 


Sl665 


Sl667 


Ratio 


Velocity 




(J2000) 


(J2000) 


(mjy) 


(mjy) 


(1667:1665) 


(kms -1 ) 


644-50 1 


53! 643 


50'.' 11 


11.4 


15.0 


1.32 


115 


1 1 -1-47 Q 


53.113 


47.88 


20.2 


<2.7 


<0.13 


47 


52.73+45.8 


52.729 


45.77 


16.9 


<2.7 


<0.16 


3 


52.42+49.2 


52.416 


49.18 


11.1 


<2.7 


<0.24 


62 


51.94+48.3 


51.937 


48.29 


11.5 


<2.7 


<0.23 


-7 


51.55+48.5 


51.543 


48.47 


<2.7 


7.42 


>2.75 


-27 


51.27+44.3 


51.276 


44.29 


13.2 


32.4 


2.44 


-113 


50.95+45.4 


50.950 


45.35 


29.1 


34.4 


1.18 


-133 


50.54+45.5 


50.542 


45.50 


5.15 


<2.7 


<0.52 


-34 


50.38+44.3 


50.381 


44.24 


22.8 


116 


5.09 


-151 


50.07+43.1 


50.066 


43.06 


<2.7 


7.66 


>2.84 


-158 


49.71+44.2 


49.714 


44.21 


<2.7 


12.1 


>4.48 


-83 


48.43+41.9 


48.432 


41.96 


6.18 


19.0 


3.07 


-125 



Table 2. Definite maser detections identified from the VLA 2006 data. Limits on the measured peak flux densities are given as 3a. 
Velocities are quoted relative to the systemic velocity of 225 km s -1 . Where a line was present in both frequencies, the measured velocities 
were the same in both frequencies. The maser IDs here are constructed from the J2000 positions measured in this dataset (relative to 
09 h 55 m +69°40') and are used throughout this paper. Maser IDs in bold are those which are new in these observations. Uncertainties 
on fluxes are ±2.7 mjy, and ±1.1 km s -1 on velocities. 
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Figure 1. All masers detected in the VLA 2006 observations at >5cr. The image is a combination of both frequencies showing masers at 
both 1665 and 1667 MHz. The contours are (—1,1,2,4,8,16,32) X 3mjy/beam and the labels correspond to the maser IDs given in Table 
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features have only one narrow peak, but a subset are broader 
with a shoulder on one side or the other. The spectra are 
discussed further in Section [5] 

Although these observations were not ideal for detecting 
absorption, some is present in the data, mainly to the east of 
the dynamical centre. The deepest absorption features are 
coincident with those found in previous Hi studies at the 
eastern end of the central kpc. 



4 NOTES ON INDIVIDUAL SOURCES 

Brief notes are given here for each maser comparing with the 
results given in lAreo et all <|2007l >. More detailed discussion 
of the properties of the population are given in the following 
section. Note that the features are identified by their J2000 
position in this dataset with labels corresponding to those 
in column one of Table [2] Where a maser was detected in 
2002, the velocities measured here are the same to within 
the uncertainties, with the exception of 51.94+48.3 (see be- 
low). Unless otherwise noted, SNRs and Hi I regions are given 



by th eir B1950 positions a s reported by McDonald et alJ 



(2002), H2O masers are from lBaudrv fc Brouilletl (Il996h and 



OlTTatelhte to features ar e from I seaouist et a! 



(Il997fl . 



Note that all masers detected in the 2002 data are also 
detected here. The only previous maser which is undetected 
here is 50.02+45.8 which was only detected in the MERLIN 
1995 observations. 

(i) 53.64+50.1 is still detected in both lines but with a 
slightly larger line ratio than in 2002. The measured ve- 
locity is the same (within the errors) for both lines and 
matches that measured in 2002 (122+9 km s _1 ). It lies 
with in a beam of the Hn region 44.93+63.9 (B1950 posi- 
tion; iMcDonald et alj|2002f ). 

(ii) 53.11+47.9 is again only present at 1665 MHz where 
its velocity (46.7+1.1 km s -1 ) is consistent with that mea- 
sured in 2002 (56±9kms _1 ). It is apparently coincident with 
the supernova remnant 44.43+61.9. 

(iii) 52.73+45.8 is also only present at 1665 MHz, as in 
2002, and the measured velocity is within the uncertainty of 
the previous result (8+9 kms -1 ). This source is co-located 
with a supe rnova remnant (44.0 1+59.6) and a satellite line 
feature (F4; ISeaauist et al1ll997f ) which is seen in absorption 
at 1720 MHz and strongly in emission at 1612 MHz with a ve- 
locity of 8+4 km s -1 , consistent with the velocity measured 
here. While listed as a supernova remnant by McDonald et 
al, the behaviour in the satellite lines led Seaquist et al to 
suggest the source as an AG N candidate, and EVN observa- 
tions by I Wills et all l| 19991 ) show comparable structure and 
brightness to SgrA. It is also within l'.'l of an X-ray poin t 
source (CXOM82 J095552.8+694047; iGriffiths et ai]|200Ch . 
but the main line rat io is atypical of megamaser sources as- 
sociated with AGNs (|Lon sdalc 200^). 

(iv) 52.42+49.2 is a new maser in this dataset. It 
is only present at 1665 MHz and is apparently coincident 
with a possible water ma ser detection (43.7+62.75) from 
iBaudrv fc Brouilletl (|l996h . listed as a possible feature only 
due to not being >5cr in more than one channel. Inter- 
estingly, the authors note that this feature has an anoma- 
lous velocity compared to the molecular emission along the 
same line of sight. The velocity measured for the OH fea- 
ture 52.42+49.2 in the 2006 data is +62.1 kms" 1 while 
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Figure 2. Spectra of each maser listed in Table[2]at 1665 MHz, arranged in order of decreasing R.A. Each maser is seen in two or more 
adjacent channels with a flux greater than five times the rms noise level. The data have been continuum subtracted and are smoothed 
with a Gaussian of width 2 channels. The a;- axis velocity scale is calculated relative to the rest frequency of the 1665 MHz line. Continued 
on next page. 
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Figure [2] continued. 
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the H2O feature of Baudry & Brouillet is at a velocity of 
— 64.6 kms -1 , putting it well off the main distribution and 
away from any other features on the eastern side of the 
galaxy on the p-v diagram (see Section \5. 411 . 

(v) 51.94+48.3 has the largest position offset compared 
to previous detections of all the masers discussed here, al- 
though the offset (0'.'36) is smaller than the size of the beam 
(1'.'4). It was the weakest detection at 1665 MHz in 2002 and 
not detected at all at 1667 MHz. In 2006 it is also only de- 
tected at 1665 MHz but with a slightly different velocity at 
it's peak (2002: 12+9 kms" 1 ; 2006: -7.3±1.1 km s" 1 ). The 
detection in 2002 was only definite at 1665 MH^] where it 
is superimposed on an absorption feature and hence was 
only just above the detection threshold. The 2006 detec- 
tion is much more significant. It is along the same line of 
sight (within the size of t he beam) as the continuum so urce 
43.21+61.3 (B1950) from iRodriguez-Rico et~ai1 (|2004T ) de- 
tected at 8.3 and 43 GHz and classified as a possible Hn 
region. 

(vi) 51.55+48.5 is a new maser in this dataset. It is 
only present at 1667 MHz where it is the weakest detection 
in this line and has no apparent associations with known 
continuum sources or other maser features. 

(vii) 51.27+44.3 is detected in both lines as it was 
in 2002. The measured position and velocity matches 
within the uncertainties (2002: -106+9 kins" 1 ; 2006: 
— 113+2 km s _1 ) but the line ratio has increased by a factor 
of two since 2002. It is within a VLA beam of both a sig- 
nificant (SI) and a possib le H2O maser (42.5+59.25) from 
iBaudrv fc Brouillet! (|l996l ). within an arcsecond of an Hn 
region (42.48+58.4) and within a beam of two others. 

(viii) 50.95+45.4 is once again also detected in both 
lines. It is the only feature where the line ratio has de- 
creased between the two epochs, but the velocities are the 
same within the uncertainties. It is apparently coincident 
with both an Hn region (42.21+59.2) and an H2O maser 
(S2). 

(ix) 50.54 + 45.5 was brighter at 1665 MHz in 2002 but 
was detected in both lines. In the 2006 observations however, 
it is only detected (weakly) at 1665 MHz. It has no known 
associations with other masers or continuum features. 

(x) 50.38 + 44.3 is apparently coincident with a known 
Hn region (41.61+57.9) and a water maser (S3). It is again 
the brightest OH maser feature detected in M82 and con- 
tinues to have the most extreme line ratio. 

(xi) 50.07 + 43.1 is a new feature in this dataset and 
is only detected at 1667 MHz. It is not coincident with any 
known continuum feature, or any other maser detection, and 
its velocity is consistent with it being located on the blue 
ward arc (see Section r6.4[) . 

(xii) 49.71 + 44.2 is again only significantly detected 
at 1667 MHz. This feature is within a beamwidth of two 
Hll regions (40.95+58.8; 40.96+57.9) and is also coincident 
with a possible H2O feature (40.9+58.25). At the same po- 
sition there is a weak feature (3.7 mjy) at 1665 MHz with a 
velocity ( — 79.8 kms -1 ) slightly offset from that measured 
at 1667MHz (—83.0 kms -1 ). However, as it is <5cr, this 
1665 MHz feature is not included in further discussions. 



(xiii) 48.43 + 1.9 is again detected in both lines and is 
brighter at 1667 MHz with a velocity comparable to that 
measured in previous datasets. It is within a beamwidth 
of an Hn region (39.68+55.6) and an H2O maser detection 
(S4). 



5 DISCUSSION 

Of the 13 features seen in these observations, ten have 
been identified in previous work i|Argo et al.ll2007l ). Five 
masers are seen in both lines, a further five are seen only 
at 1665 MHz and three are present only in the 1667 MHz 
data. The only maser reported previously which has not 
been detected again here is 50.02+45.8, a source which was 
only detected in low-velocity resolution MERLIN observa- 
tions carried out in 1995 and, although it was detected at 
>3a in 1995, it was the weakest detection in that dataset 
and has not been detected subsequently. 

This section discusses apparent associations, spatial 
structure, line ratios, spectra, velocity distribution and the 
blue ward arc in more detail. 

5.1 Possible associations 

Of the five masers visible only at 1665 MHz, two are ap- 
parently coincident with supernova remnants (53.11+47.9; 
52.73+45.8), one of which is also co-located with a satel- 
lite line feature (52.73+45.8), one lies within a beam of a 
possible H2O feature (52.42+49.2), and two have no ap- 
parent associations (51.94+48.3; 50.54+45.5). Of the three 
sources seen only at 1667 MHz, two are new detections 
(51.55+48.5; 50.07+43.1) neither of which appear associ- 
ated with known continuum features or other maser detec- 
tions. The third of these, 49.71+44.2, was previously noted 
to coincide with an Hll region. The five features visible in 
both lines (53.64+50.1; 51.27+44.3; 50.95+45.4; 50.38+44.3; 
48.43+41.9) are all apparently coincident with known Hll 
regions. 

As has been noted previously, considering that the disk 
of M82 is viewed almost edge on these apparent associations 
are likely to be due to line of sight effects. These observations 
lack the required spatial resolution to determine whether the 
sources are physically related. 

5.2 Spatial structure 



1 Figure 4 of lArgo et al. (2007) shows the incorrect spectrum for 
this maser (labelled 51.87+48.3). 
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Figure 3. Spectra of each maser listed in Table[2]at 1667 MHz, arranged in order of decreasing R.A. Each maser is seen in two or more 
adjacent channels with a flux greater than five times the rms noise level. The data have been continuum subtracted and are smoothed 
with a Gaussian of width 2 channels. The x-axis velocity scale is calculated relative to the rest frequency of the 1667 MHz line. Continued 
on next page. 
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Figure [3] continued. 



As with the 2002 VLA data, the resolution of these 
observations is insufficient to see spatial structure in the 
maser regions. This is not surprising given that the physical 
size of typical Galactic masers placed at the distance of M82 
would be much smaller than the linear resolution of the VLA 
at that distance (22 pc at f.6GHz). 

The only maser which shows convincing physical exten- 
sion is 50.95+45.4 near the centre of the galaxy. Previous 
VLA observations also showed a similar extension to the 
SW which was best fitted by three Gaussians. Higher res- 
olution observations with MERLIN sho wed a weak secon d 
component along the same position angle l|Argo et al.ll2007l ). 
These observations show an extension in the same direction 
(Fig. [1]) which varies in frequency: the position of the peak 
moves from NE to SW as frequency increases. The same ef- 
fect can be seen in both main lines and can also be seen 
clearly in the p-v diagram (see Section \5A\ . It is likely that 
the extension of 50.95+45.4 is due to several maser spots 
lying fairly close together but with different velocities, the 
resolution of the VLA at this frequency is not sufficient to 
distinguish them. 

In the 1667 MHz data, 51.27+44.3 appears to be slightly 
extended to the east. Again, due to the small size of known 
Galactic maser regions, this is likely caused by the superpo- 
sition in space of two unresolved maser spots, rather than 
any intrinsic structure. This is also illustrated by features 
50.54+45.5 and 50.38+44.3 which are located very close to- 
gether on the sky, but are clearly separate features as their 
velocities are —34 and — 151 kms -1 respectively. 



5.3 Line ratios 

For an optically thin cloud in local thermal equilibrium, the 
expected ratio for emission from the 1665 and 1667 MHz 
lines is 1.8. Column 7 of Table [2] shows line ratios of each 
feature present in this dataset. Where a feature is only de- 
tected in one line, the limit is calculated from the 3<r noise 
in the other line. Note that, compared to typical velocity 
resolutions used in observations of Galactic masers, these 



observations are still under-resolving the masers in M82, so 
it is possible that they could be brighter still. 

In 2002, six features were detected in both lines al- 
lowing line ratios to be calculated. One of these features, 
50.54+45.5, had a line ratio of 0.69 in 2002 but is not de- 
tected above 5<r at 1667 MHz in 2006. The other five masers 
in this group were all detected above 5<r in more than two 
channels in the 2006 observations in both lines. Of these, 
three have line ratios which are consistent with those mea- 
sured in 2002, given a typical 5 per cent flux uncertainty 
on each measurement. However, 51.27+44.3 and 48.43+41.9 
both have ratios significantly larger than those measured in 
the 2002 data. In the case of 51.27+44.3 the ratio in 2006 is 
twice that measured in 2002. 

As in 2002, 49.71+44.2 is again only detected at 
1667 MHz. Two of the three new maser features (51.55+48.5 
and 50.07+43.1) are also only present in the 1667 MHz 
dataset and so have limits calculated from the 3a noise 
measured in the 1665 MHz line. 53.11+47.9, 52.73+45.8 and 
51.94+48.3 are again only seen at 1665 MHz in this dataset. 
50.54+45.5 was previously detected at >5<r in both lines in 
2002 where it was stronger at 1665 MHz, but in the 2006 
data it is only visible at 1665 MHz. 

With the exception of 50.95+45.4, for all masers which 
are detected in both lines the ratios have increased since the 
2002 observations. The greatest difference is for 51.27+44.3 
which has doubled from a ratio of 1.2 in 2002 to 2.4 in 2006. 
The feature 50.38+44.3 continues to have the most extreme 
line ratio, 4.7 in 2002, 5.1 in 2006. Interestingly, while it 
is a significant distance from the dynamical centre of M82 
and apparently coincident with an Hn region, 50.38+44.3 
has a high line ratio more typical of megamasers associated 
with AGN than Galactic masers in regions of star formation 
jLonsdald^Ool l. 

One caveat with comparing measured line ratios be- 
tween data sets is that whether or not the line is resolved 
affects the measured peak flux. In some cases, the lines are 
clearly resolved in frequency, whereas with previous mea- 
surements they were unresolved. 
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Figure 4. Position-velocity plots for each IF in the 1665 MHz dataset. Features are labelled according to their positions as listed in 
Table[2] Contours are (—10 to 10) X 1.6mjy/bm. The velocity Okms" 1 represents the systemic velocity of M82 (225 kms ) and the 
centre of M82 is taken to be 09 h 55 m 52?132 +69°40'46'/14 (J2000). 




Figure 5. Position-velocity plots for each IF in the 1667 MHz dataset. Features are labelled according to their positions as listed in 
TablefS] Contours are (-10 to 10) X 1.6mjy/bm (left) and (-10 to 10, then 12, 16, 24, 40, 68) X 1.6mjy/bm (right). 
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5.4 Spectra and velocity structure 

Figures [2] and show the spectra of each feature >5<j in 
more than one consecutive channel at 1665 and 1667 MHz 
respectively. The spectra are arranged in order of decreas- 
ing R.A. (left to right in Figure [TJ and labelled according to 
their IDs in Table [2] The p-v plots in Figures [4] and [5] show 
the distribution of masers in velocity along the major a xis of 
the galaxy. Following the method of IWills et al.l (|2000T ). be- 
fore producing the p-v plots each cube was first rotated 17° 
anticlockwise so that the major axis of the galaxy was hor- 
izontal. Then each cube was compressed to one plane along 
the declination axis, ignoring the five channels at either end 
of the band where the response was poor. 

The velocity distribution of the masers in this dataset 
lies along the same di stribution as that seen in the CO 
and Hi distributions in IWills et all (|200Qh . As seen in the 
2002 observations however, several masers lie along an arc 
blue ward of the main distribution. As with the Hi gas, the 
strongest features of this arc lie on the same axis as that of 
the io nised [Ne n] gas distribution of lAchtermann fc Lacvl 
(| 19951 ') . The cause of this feature has been suggested to be 
eith er an expanding su perbubble or the xi orbits of an inner 
bar (|Wills et alJl2000D . 

The spectra of the two brightest masers, 50.95+45.4 
and 50.38.44.3, both show a blue wing at 1667 MHz which 
is less obvious at 1665 MHz. It can be seen from the right- 
hand panels of Figures [4] and [5] that this region of the p-v 
diagram has some interesting structure. Both 50.95+45.4 
and 50.38+44.3 show velocity structure with significant ex- 
tension in both lines, structure which is more extended at 
1667 MHz than it is at 1665 MHz. This can be seen more 
clearly in Figure [6] where the two p-v plots are superim- 
posed for this region. The redward extension of 51.27+44.3 
and the low-level bluest third peak of 50.38+44.3 are only 
visible at 1667 MHz, while 50.95+45.4 has broadly the same 
structure in both lines with extensions to both the east (red) 
and west (blue) and a double-peaked structure, with the ex- 
tension again greatest at 1667 MHz. 

Another maser which shows velocity structure is 
50.07+43.1 which appears to have two narrow peaks, but 
is only present at 1667 MHz. The two peaks here are more 
closely spaced in frequency than those of 51.27+44.3 and 
have similar brightnesses, with the bluer peak marginally 
stronger. 

In the p-v diagram, the feature 51.27+44.3 also clearly 
has two peaks at 1667 MHz with the bluer peak significantly 
stronger. The p-v plots in Figures [4] and [5] show that there 
is also emission at the same position on the major axis but 
with an even bluer velocity, blending into the large fea- 
ture at 50.95+45.4. This illustrates the complexity of the 
ISM in M82. The second, weaker peak in the spectrum of 
51.27+44.3 in Figures [2] and [3] is actually this weak fea- 
ture close to 50.95+45.4. Given the linear resolution of the 
VLA, it is likely that this feature is a different maser spot, 
physically unconnected but along the same line of sight as 
51.27+44.3, with a velocity similar to that of 50.95+45.4. 

The structure of this group of masers in the p-v diagram 
suggests several closely spaced maser regions sited along the 
blue ward arc, and it is likely that this region contains many 
clouds at different velocities spread out over a significant 
volume. This region is discussed further in Section \6. 41 
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Figure 6. The region around 50.95+45.4 in more detail. The 
1665 MHz data is shown as greyscale with a range from to 
25 mJy/bm, while the 1667 MHz data is plotted as contours of 
(-1,1,2,4,8,16,32) X 1.9 mJy/bm. The emission at 1667 MHz is no- 
ticeably broader in frequency than that at 1665 MHz. 



The only source located in the main velocity distri- 
bution to show significant velocity structure is 48.43+41.9 
which is weak but has a double-peaked profile in the p-v plot 
at 1665 MHz, but a stronger single peak with pronounced 
wings at 1667 MHz. 



6 COMPARISON WITH OTHER 
WAVELENGTHS 

Due to its proximity, M82 has been well studied across the 
electromagnetic spectrum. It is interesting to compare the 
emission seen here in OH with features seen at other wave- 
lengths in order to gain a greater understanding of the con- 
ditions and dynamics within the central kiloparsec of the 
galaxy. 

6.1 Stars and dust 

Figure [7| shows the positions of the masers plotted on an 
archival image from the ACS instrument on the HST using 
a 680-980 nm filter. While it should be noted that there are a 
lack of directly comparable features between the optical and 
radio fields, the astrometry of the HST image was tweaked 
using positions of objects in the field from the 2MASS and 
SDSS catalogues, with a resulting estimated rms positional 
uncertainty of 0'.'6 (the symbols are 1 arcsecond in diam- 
eter). Different symbols indicate whether the maser is de- 
tected at just 1665 MHz, just 1667 MHz, or both. As this 
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Figure 7. The OH masers in M82 plotted over the optical emission from the central starburst. Grey scale: Ho HST ACS image (program 
ID 9788, PI Ho). The positions of the masers are also plotted. Those visible at just 1665 MHz are shown as plusses, those detected only 
at 1667 MHz as crosses, and those present in both lines as stars. 



figure shows, and as might be expected from maser emission 
associated with star formation regions, most of the masers 
are located within dusty, heavily obscured regions of the cen- 
tral starburst, particularly the masers located on the blue 
ward arc. 



6.2 Atomic and molecular gas 

To carry out a comparison with the absorption seen in the 
atomic gas traced by Hi absorption, 1420 MHz data from the 
VLA archive has been used from an A-ar ray observation o n 
November 27th 1996 (programme AW444; lwTlls et al.ll2000h . 
The dataset was reduced using standard spectral line proce- 
dures in AIPS and a p-v plot was produced using the same 
method as described in Section[5]4] While the spatial resolu- 
tion of the two datasets is similar, the 1996 observation was 
carried out with a bandwidth of 3.125 MHz over 63 channels 
resulting in a velocity resolution of 10.3 km s -1 . To com- 
pare this with the higher velocity resolution OH data from 
2006, the two IFs at each frequency in the OH observations 
were combined using the AIPS task UJOIN resulting in two 
datasets of 190 channels each, one for each of the OH lines. 



The Hi data has only 63 channels so to compare the datasets 
the number of channels was artificially increased using the 
AIPS task xsmth which can interpolate over channels. The 
resulting p-v plots are shown in Figure [8] 

The deep absorption in the OH data is coincident with 
the deepest absorption features seen in the Hi data at the 
eastern end of the disk. The masers appear to avoid the 
deepest Hi absorption features, but this could be due to OH 
absorption being present at the same locations since strong 
absorption can hide weak maser peaks (as illustrated by the 
2002 VLA data). 

The maser features on the blue arc are coincident with 
the arc seen in the Hi data on the blue war d side of the 'hole ' 
ilWills et al.ll2000l) also seen in the C O i|Shen fc Lol Il995l ) 
and [Ne n] (jAchtermann fc Lacvlll995l ). This feature is dis- 
cussed further in Section 16.41 More recently, ISeaauist et al.l 
(2006) mapped M82 in 12 CO J=6-5 and found evidence for 
the same blue feature, and a corresponding weak feature to 
the east of the galaxy, as well as higher than average excita- 
tion in the region surrounding 41.95+57.5 where the bright 
masers on this feature are located. 
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Figure 8. A comparison between the OH and Hi datasets. The 1665 MHz p — v plot is on the left, 1667 MHz on the right. The greyscale 
shows the Hi absorption while the contours show the OH absorption and emission. The greyscale range is from to —28 mjy/bm, 
while contours are drawn at (-4,-2,-1,1,2,4,8,16,32,64) X 2.4mjy/bm (left) and 1.5mjy/bm (right). The centre of M82 is taken to be 
09 h 55m52f 132 +69°40'46'.'14 (J2000) and the velocity has been corrected for a systemic velocity of 225 kms -1 . 



6.3 Ionised gas 

Comparin g the positions of the masers with the [Nell] emis- 
sion from lAchtermann fc Lacvl (1 19951 ). several masers ap- 
pear around the edge of the 12.8/im ring. The three bright 
features around the location of the radio continuum source 
41.95+57.5 also appear to coincide with the [Nell] peak la- 
belled Wl. This is not surprising given that their velocity 
distribution also coincides with that of the [Nell], rather 
than with the main disk rotation; the velocity distribution 
of the [Nell] emission is much steeper than that of the main 
disk traced by the Hi or OH absorption and the western 
side of this steeper distribution is well traced by the brighter 
masers in this dataset. 

Figure [9] shows the positions of the masers relative to 
the ionised gas as traced by 408 MHz radio emissio n from 
a MERLIN observation in 1994 (I Wills et al.lll997h . From 
this figure it can be seen that the masers tend to avoid the 
peaks of low frequency emission and, in particular, several 
of the masers are located within the 'hole' feature noted 
by Wills et al. This feature is ~100pc in diameter, roughly 
centered on the compact object 41.95+57.5, and has been 
suggested to be due to free-free absorption by gas in a gi- 
ant Hll region ionised by a large stellar cluster (note that 
this is also the locat ion of the high excit ation seen in 12 CO 
J=6-5 transition bv lSeaquist et al.ll2006l mentioned above). 
The masers located in this region however are not typical 
of Galactic masers found in star forming regions as they all 
(except 50.54+45.5) have line ratios >1. Interestingly, it is 
also this group which are located on the blue feature in the 



p-v diagram except (again) 50.54+45.5 which has a velocity 
which places on the main distribution. 



6.4 The blue arc 

As has already been noted, an obvious feature of the p-v 
plots is the arc to the west of the disk which lies blue ward 
of the main distribution. It has been shown previously that 
the OH absorption traces the arc well, and that several of 
the brightest mas e r features in M82 are located on the arc 
jArgo et alj|2007t ). lwills et al.l ( |2000h compare the distribu- 
tion of Hi traced by absorption with that of CO emission and 
ionised [Nell] finding a similar feature in both Hi and CO, 
and noted that the steeper gradient of the [Nell] gas matches 
at least the eastern edge of this feature in the Hi and CO 
distributions. This informa tion was used to model a bar fea- 
ture in the centre of M82. iMatsushita et ail (|2005l ) showed 
a p-v plot of 12 CO(1-0) together with H41a and noted that 
the known OH and H2O masers were mainly clustered in 
velocity around the "superbubble" . This arc is also present 
in recent op tical data of emission lines over the central kpc 
of the disk l|Westmoquette et alJkoogl ). 

The gradient of the main distribution measured from 
these data is found to be ~6kms _1 arcsec , in agree- 
ment with ~6.5k ms~ 1 from Hi and CO measurements of 
IWills et alj (|2000T ). whereas the gradient of the emission on 
the eastern edge of the blue ward arc is ~20 km s _1 arcsec - , 
similar to that of the [Nell] gas distribution. The three 
brightest masers at 1667 MHz are all located on this feature 
while the fainter features are on the main distribution ex- 
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Figure 9. Positions of the masers plotted over the 408 MHz emission from MERLIN observations in 1994 from IWills et al ] lll997ly 
Contours are plotted at (-2, 2, 4, 6, 8, 10, 12, 14, 16, 18, 20) X lmjy/bm. As for Figure [JJ masers visible at just 1665 MHz, just 
1667 MHz, or both, are shown by plus, cross and star symbols repectively. The masers avoid the areas of emission with the brightest 
masers sited in the extended absorption region surrounding the source 41.95+57.5, see Section 16.31 



pected for a rotating disk. With the exception of 50.54+45.5, 
all of these masers are weaker at 1665 MHz. 

50.54+45.5 is the odd one out in this group since, de- 
spite appearing to be part of the group spatially, dynam- 
ically it is quite separate. Velocity measurements clearly 
place it not in the blue arc on the p-v diagram but on 
the main disk distribution tracing the bulk rotation of the 
galaxy. It is also the only maser in this group which is visi- 
ble only at 1665 MHz, giving it a line ratio of <1 and, of the 
masers in this group, 50.54+45.5 lies closest to the edge of 

the dust lane in the optical map. 

These observations confirm the results from I Argo et al.l 

that the brightest masers exist on this blue ward fea- 
ture. As well as being clustered in frequency, they are clus- 
tered spatially, in an apparent arc to the north of the strong 
continuum source 41.95+57.5. Comparing the positions of 
the masers with optical emission in Figure [7] the masers 
in this group are located in a thick dust lane to the South 
of a bright knot. They are also all spa tially located around 
the [N ell] features Wl and W2 from lAchtermann fc Lacvl 
|l995l ) which shows a similar velocity distribution to that of 
the masers. 



7 CONCLUSIONS 

Sensitive, high spectral resolution observations of the M82 
starburst have been carried out at 1.6 GHz using the VLA 
in A-configuration. The observations confirm previous de- 
tections and discovered three previously undetected masers. 
All the masers are significantly brighter than typical Galac- 
tic masers, an effect which is probably due to the super- 
position of several masing clouds along a particular line of 
sight, a conclusion supported by the velocity structure we 
are starting to see in these higher resolution observations. 



The spatial distribution of the masers follows that of the 
continuum emission with most of the masers apparently co- 
incident with continuum objects and/or other maser species. 
Their distribution in p-v space largely follows that of the ab- 
sorption seen in Hi. Most of the masers are located on the 
main velocity distribution which traces the bulk motion in 
the disk, however three of the brightest masers are located 
on a significant blue ward feature around a 'hole' observed 
previously in Hi and CO, and traced particularly well by 
[Nell] emission. If this group of masers is tracing out an in- 
ner ring, as suggested by comparisons with the [Nell] emis- 
sion, then it might be expected that a corresponding group 
would be seen on the eastern side of the dynamical centre 
with redshifted velocities compared to the bulk distribution. 
Such a feature is not traced either by maser emission or OH 
absor ption, although a h int of such a structure was noted in 
Hi bv lWills et al] feooch . 

Despite the improvement over previous observations, 
many of the masers are still largely unresolved in frequency. 
However, those on the blue ward feature do show some in- 
teresting structures with multiple velocity components, sup- 
porting the idea that several maser features are being de- 
tected along one line of sight. This structure is apparent 
at both 1665 and 1667 MHz with all three masers more ex- 
tended on the p-v diagram at 1667 MHz. 50.95+45.4 shows 
particularly interesting structure with evidence of multiple 
components both spatially and in velocity. 
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